Inter-dependence of the volume and stress ensembles and equipartition 
in statistical mechanics of granular systems 
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We discuss the statistical mechanics of granular matter and derive several significant results. First, 
we show that, contrary to common belief, the volume and stress ensembles are inter-dependent, 
necessitating the use of both. We use the combined ensemble to calculate explicitly expectation val- 
ues of structural and stress-related quantities for two-dimensional systems. We thence demonstrate 
that structural properties may depend on the angoricity tensor and that stress-based quantities 
may depend on the compactivity. This calls into question previous statistical mechanical analyses 
of static granular systems and related derivations of expectation values. Second, we establish the 
existence of an intriguing equipartition principle - the total volume is shared equally amongst both 
structural and stress-related degrees of freedom. Third, we derive an expression for the compactivity 
that makes it possible to quantify it from macroscopic measurements. 

PACS numbers: 64.30. -|-t, 45.70.-n 45.70.Cc 
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The statistical mechanical formalism, introduced to de- 
scribe granular materials [1-3], was expected to be a plat- 
form for derivations of experimentally measurable equa- 
tions of state and constitutive relations. It has not yet 
lived up to its full potential due to several difficulties 
that traditional thermodynamic theories do not suffer 
from: lack of ergodicity, uncertainty over the identities 
and number of degrees of freedom (DoF), and the dif- 
ficulty to realise a simple analog of a thermometer - a 
'compactometer'. Whilst these problems reflect more on 
the application of the theory rather than on the theory 
itself, a more serious concern has arisen from recent sug- 
gestions of an absence of an equipartition principle[4, 5] 
in agitated systems. Here we derive a number of signifi- 
cant results. First, we show that the correct phase space 
consists of both structural and force DoF, thus calling 
into question much of the results in the literature, ob- 
tained from either ensemble alone. Second, we show the 
existence of an equipartition principle in two-dimensional 
static systems. Third, we show that, in such systems, the 
compactivity can be quantified directly from macroscopic 
mean volume measurements. 

The initial statistical mechanical approach was based on 
a volume partition function of N{» 1) grains[l]. 
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where W is a volume function that sums over all the 
possible volumes that basic volume elements can realise 
and Xq is the compactivity - a measure of the fluctua- 
tions in the ensemble of realisations that is the analogue 
of temperature. The structural DoF (SDF), identified 
explicitly below, are all the independent variables that 
determine the structure of an assembly of grains in me- 
chanical equilibrium, given the mean number of force- 
carrying contacts per grain z[6]. 



The partition function Zv enables almost closed ther- 
modynamics. For example, once W and the SDF have 
been identified (see below), the mean volume 2D can 
be computed. Nevertheless, Zy is unable to specify the 
macroscopic state of the system completely, since the en- 
tropy remains only partially accounted for, it leaves out 
an entire set of microstates - those of different stress 
states. These microstates are described by a differ- 
ent partition function, Z/[7-9], an idea supported later 
numerically [9, 10]. The stress ensemble gives rise to the 
partition function 



Zf 



'djall boundary forces} (2) 



Here a, (3 run over the Cartesian components x, y and 
J^ajS are the components of the force moment function, 
from which the stress CTa/j is derived. 
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Here the sum runs over pairs of grains in contact gg', 
TZP^ is the position of the contact point, measured from 
the centroid of grain g, F^^ is the force that g' applies to 
g and is the volume associated with grain g. = 
dJ'ap/dS has been named 'angoricity' - a tcnsorial ana- 
logue of the temperature and the compactivity [7, 11, 12] 
- and S is the entropy, defined as the log of the number 
of both structural configurations and stress states. Note 
that integrating this partition function, any expectations 
values would be a function of all TZ^^ . 
The volume and stress ensembles are treated in the lit- 
erature as independent, leading to a total partition func- 
tion Z = ZyZj. Consequently, results have been derived 
solely from the statistics of one ensemble or the other. 
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Here, we challenge this view. We argue that such deriva- 
tions are misguided and we outline the calculation of the 
correct partition function. Using the correct ensemble 
we demonstrate derivation of a number of expectation 
values in two dimensions (2D), including the expected 
intergranular force distribution, and derive a surprising 
equipartition principle for static systems. 




We put forward the following three arguments. 

1. The volume ensemble alone is insufficient to describe 
the entropy of mechanically stable granular systems. 

A volume ensemble implies the exact same boundary 
forces. However, many-grain experiments cannot repro- 
duce the same grain configurations, nor the precise forces 
on every boundary grain. Only global boundary stresses 
can be controlled, i.e. averages over boundary force com- 
ponents. Thus, the statistics of the boundary forces must 
be taken into consideration. 

2. The stress ensemble alone is insufficient to describe 
the entropy of mechanically stable granular systems. 
The stress ensemble comprises a fixed granular configu- 
ration, to which all combinations of boundary forces are 
applied. The ensemble is subject to constraints, e.g. that 
the total boundary stresses are fixed. Such a system can- 
not be realised experimentally in very large assemblies 
(albeit possible in numerical experiments). Indeed, any 
integration over Z f remains a function of the SDF. 

3. The volume and stress partition functions are inter- 
dependent, Z ^ ZvZf. 

This statement follows from the above two arguments - 
correct calculations of expectation values must be based 
on a combined ensemble of all structural arrangements 
and all boundary forces. Specifically, this is a conse- 
quence of the explicit dependence of both the volume 
function in Z^ and the force moment function in Z f on 
the structural DoF (SDF). 

The above arguments hold in any dimension and we pro- 
ceed to illustrate them explicitly in 2D. Consider an en- 
semble of 2D iV-grain systems (A^ ^ 1), each of mean 
contact number z. The systems are in mechanical equi- 
librium under M external compressive forces, acting on 
the boundary grains. We disregard body forces, in the 
absence of which 'rattlers' can also be ignored, as they 
do not affect the stress states in static piles. 
It was proposed to use 'quadrons' [13, 14] as the ele- 
mentary volumes, both in two and in three dimensions. 
These are space-tessellating (generically) quadrilateral el- 
ements (figure 1). The quadron is constructed on two vec- 
tors as its diagonals: connects contact points around 
a grain in the clockwise direction and R'^ extends from 
the centroid (mean position) of the contacts around the 
grain to the centroid of the contacts around a neigh- 
bour cells. In terms of these, the volume function is 
W = '}2,qV'^ = ^|r^ X i?'! (summation implied over re- 
peated indices) and the partition function is 



FIG. 1. Tlie vectors r^ connect contact points clockwise 
around grain g. The vectors R' connect from grain g cen- 
troids to cell c centroid. These vectors are the diagonals of 
quadron q (blue). 
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Here Nsdf is the number of SDF, discussed below. The 
vectors can be expressed as linear combinations of 
the [8] and, since the latter close loops, only Nz/2 of 
them are independent [6, 8], leading to Ngdf = Nz. Defin- 



ing the vector p = (^rl , 
becomes exactly quadratic and we have 
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(5) 

Here p, q run over quadrons, a, /3 run over vector compo- 
nents X, y and yl is a matrix whose elements are 



xy 

"'yy 



q,p < Nz/2 
q < Nz/2 , p > Nz/2 
p < Nz/2 , q > Nz/2 
q,p > Nz/2 



We can now evaluate Zy. Assuming a uniform measure 
of the DoF and that the contribution of very large r mag- 
nitudes is negligible, (5) can be calculated explicitly 



Zv 



\2n) 



Nz 



detA 
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The stress partition function consists of all the possi- 
ble combinations of compressive forces on the bound- 
ary grains, (f" (to = 1,2,...,M), subject to the con- 
straint that the total stress on the boundary is fixed 
[6, 11, 12]. Since the configuration is presumed fixed, 
only boundary forces that do not drive the system out of 
mechanical equilibrium are permissible, i.e. the bound- 
ary stresses must be below the yield surface. It has been 
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FIG. 2. fa is the a component of the loop force f^, which 
contains the quadron q (shaded blue). The loop forces of c 
and c' give the inter granular force at the contact point that 
they share, F®® — f ~ f [13]. Quadron p (striped) shares 
the same loop as q and hence also the same loop force. 



shown[13, 14] that the force moment function (3) can be 
written as 



(7) 



where is the a component of a loop force of the 
cell containing the quadron q (see figure 2). The loop 
forces are defined in terms of the contact forces[13], e.g. 
pgg = — figure 2, and they conveniently satisfy 

force balance conditions. In 2D assemblies of arbitrarily- 
shaped frictional grains, there are N (z/2 — 1) loop forces 
(a consequence of Eulcr relation), of which N/2 can be 
determined from the torque balance conditions. For clar- 
ity, we next focus on isostatic systems (z = 3) ; extension 
to hyperstatic assemblies (z > 3) is possible. Substitut- 
ing (7) into (2) gives 



Zf = I e 



-<j}-E^ B p 



M 

lid 



M 



g-C' -E^ B p 



(10) 

To compute the total partition function, we need to in- 
tegrate over the combined volume-stress phase space, 
dZ = dZydZf 



\p-A-p-g-Q' 



{d^'-p) {d'^'g) (11) 



where we have defined, for brevity, Q =B^-E-C. This 
expression establishes our claim that Z is not the product 
of Zv and Zf of eqs. (6) and (10). 

Integrating (11) is straightforward due to the integrand's 
Gaussian form and we use it next to calculate several 
expectation values. The exponential contains a linear 
and a quadratic term in p and, completing to quadrature 
and changing variables to p = p+A^^Qg, we can separate 
the variables in the exponent, 



where we define for shorter notation P = Q'^ ■ A^^ ■ Q. 
Calculating the mean volume we obtain: 



Xn 



[pAp + gPg)e^^ 



which separates into two gaussian integrals, giving 



M 



II dr 



(8) 



where the integration runs over all the independent 
boundary forces g^. Since quadrons sharing the same 
cell have the same loop force (figure 2) and the loop 
forces depend linearly on the M boundary forces, then 
only N/2 of the Nz quadron forces are independent. It 
is therefore convenient to define the loop forces vector 

= (/I , /J , . . . , /f /^ /I , /2 , . . . , A^/') . The solution for 
(j) in terms of the boundary forces can be written as 



1/3 9/3 



(9) 



where a, ^ = x,y, c = l,2,...,iV/2 runs over all cehs, 
TO = 1, 2, . . . , M runs over all boundary forces and C is 
N X 2M matrix. In terms of these, = Ecj), where E is 
a Nz X N matrix. Defining further i?^^ 



6qp being the delta function, we finally obtain Zf 



(13) 



This result is significant for several reasons. First, it 
is independent of the details of the connectivity matrix 
A and of the particular stress state. Second, it reveals 
a striking equipartition principle: the mean volume is 
shared equally among the zN structural and the 2M 
force DoF, each having on average Xq/2. It is analo- 
gous to the mean energy per DoF in thermal systems 
3fcBT/2, but we emphasize that no energy is involved 
in this formalism. Third, it quantifies the compactiv- 
ity in terms of measurable quantities. An important 
consequence of this finding is that it makes possible to 
start analyses 'inductively' by assuming that the volume 
per DoF is Xq/2, as done as standard in thermal sys- 
tems analogously with /cbT. Note that using only Z„ 
(eq. 5) gives {Vy) = zNXo/2, which overestimates the 
compactivity. 

All relevant expectation values can be expressed in terms 
of p and g and hence evaluated, albeit with more algebra. 
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2M 

E 



dR,. 



2M 

E 



(14) 
(15) 
(16) 



where R ^ ■ ■ ^ ■ Q ■ Y, T = ■ ■ 
A-'^ ■ A-^ ■ Q -Y , U = Y^ ■ ■ E'^ ■ E ■ C - Y ,Y ]s the 
matrix that diagonahses P, pi are the eigenvalues of P, 
and are straightforward functions of E and C . Note 
that both resuhs (13) and (16) are directly relevant to 
experimental measurements [15, 16]. These exact results 
do more than demonstrate the utility of the combined 
ensemble, they reveal unexpected dependences of these 
expectation values on the compactivity and angoricity. 
For example, {p ■ p) is not only proportional to Xq, as 
one would expect, but it also depends on X^p via a ho- 
mogeneous function (HF) of order 0. Also unexpectedly, 
the mean inter-granular force magnitude square (/ • /) is 
both a HF of order 2 of Xa/j and linear in Xq. Yet {J^ap) 
is, unsurprisingly, a HF of order 1 of the angoricity and 
independent of Xq. These results show the significance 
of using both the stress and the volume ensembles. 
To conclude, we have presented three main results. First, 
we have shown that the compactivity-based volume en- 
semble and the angoricity-based stress ensemble are de- 
pendent and need to be used simultaneously. We reiter- 
ate, the entropy is the log of all the microstates, which 
include both the SDF and the stress states - because Zy 
and Zf are dependent, it is not simply the sum of the con- 
figurational and stress entropies. This calls into question 
the large body of work, obtained from cither ensemble 
alone. We have used the combined partition function to 
obtain explicitly the expectation values of: the mean vol- 
ume, force moment, distance between intra-grain neigh- 
bour contact points, and contact force magnitude. We 
find, surprisingly, that (y) depends explicitly on the force 
degrees of freedom, that (p • p) depends on the angoricity, 
and that (/ • /) on the compactivity. Second, the calcu- 
lation of (y) reveals the existence of an equipartition 
principle - the mean volume of static systems is shared 
equally amongst both structural and force-related DoF, 
with each getting a volume of Xo/2. This result shows 
that, although equipartition is questionable in dynamic 
dense systems[4, 5], it exists for static ones. Moreover, 
since static granular systems are the equivalent of "zero 
temperature" granular fluids, this result gives hope that 
an equipartition principle may be found for dense dy- 
namic systems by extending dynamic descriptions to in- 
clude structural and force DoF. Third, we have derived 
an expression for the compactivity in terms of measur- 



able quantities - the mean volume and the mean contact 
number and the loading forces. A significant implication 
of our results is that the compactivity and angoricity are 
not in themselves the conjugate variables of volume and 
force moment, as previously believed. Instead, it is the 
expression in eq. (12) that represents a convolution of the 
volume and force moment functions with the compactiv- 
ity and angoricity. 

It would be interesting to test our analysis experimen- 
tally and numerically, e.g. by constructing assemblies 
at different compactivities and angoricities and examin- 
ing expectation values as functions of these parameters. 
Furthermore, since the arguments establishing the inter- 
dependence of the volume and stress ensembles hold in 
any dimension, it should be possible to extend our anal- 
ysis to 3D systems, at least numerically. 
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